Maxwell stress in fluid mixtures 
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We examine the structure of Maxwell stress in binary fluid mixtures under an external electric 
field and discuss its consequence. In particular, we show that, in immiscible blends, it is intimately 
related to the statistics of domain structure. This leads to a compact formula, which may be useful 
in the investigation of electro-rheological effects in such systems. The stress tensor calculated in a 
phase separated fluid under a steady electric field is in a good agreement with recent experiments. 
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Introduction — Like any other complex fluid, blends of 
immiscible fluids exhibit rich phase behaviors and dy- 
namics under external fields [H. Their domain interface 
is inherently soft so that the configuration is easily de- 
formed, from which a large stress contribution shows up. 
Rheological properties under flow fields have been well 
understood at a semi-quantitative level [2| where the so- 
called interfacial tensor plays a central role to capture 
the statistical properties of complicated interconnected 
domain structures. 

The shape change of droplets and interfacial instabili- 
ties under an electric field in immiscible fiuids having dif- 
ferent dielectric constants have been studied both exper- 
imentally and theoretically for many years Struc- 
tural transitions induced by electric field have also been 
investigated in microphase separation in block copoly- 
mers |7H9|. However, these previous studies were con- 
cerned mainly with the morphological change of domains. 
It should be emphasized that a change of domain struc- 
tures drastically affects flow behavior of the system and, 
therefore, produces a unique rheological effect. This type 
of electro-rheology is of fundamental importance since 
the cross coupling between flow fleld and electric fleld 
is relevant to characterizing the departure from equilib- 
rium. Recently, experiments of electro-rheology in phase 
separating fluids have been conducted A key 

physical quantity is the electric (Maxwell) stress, which 
is expected to be intimately related with the spatial do- 
main structures. 

In this paper, we discuss a fundamental relation be- 
tween the Maxwell stress and domain conflgurations in 
fluid mixtures. We start with the basic equations for 
phase separation dynamics in which the free energy 
functional contains the electro-static energy. First, we 
show that the reversible mode coupling term in the dy- 
namic equation for the local velocity produces exactly 
the Maxwell stress. We then eliminate the local elec- 
tric fleld E{r) and obtain an expression to the Maxwell 
stress in terms of the external electric field E"^'^ , which 



involves a non-local coupling among concentration fiuctu- 
ations of the induced dipole type. Then, for immiscible 
blends, the use of the basic statistical property of ran- 
dom configuration of interfaces at short distance (known 
as the Porod law) makes us propose that the macroscopic 
Maxwell stress contribution can be represented in terms 
of the interfacial tensor. This provides us with a useful 
formula which connects the Maxwell stress with arbitrary 
configuration of domain structures. 

Basic equations — We consider a binary fiuid of A and 
B components, the local volume fractions of which are 
represented by <^yi(r) and (j)B{r). The free energy for 
(j}{r) = 4>A{r) — (psir) consists of two parts: 



F{^ir}} = i^i{0(r)} + F2{0(f), E{r)}, 



where 



Fi{H^} = / dr 



(1) 



(2) 



and 



J^2{0(r),i?} 



1 
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dr e[<j>{r)]E{r) ■ E{r) (3) 



The constant K is positive, /(</)) is a polynomial of cj) 
with two minima. The second free energy F2 is the elec- 
tric energy where E{f^ is an electric field. The dielectric 
constant e is assumed to depend on as e = e -I- 5e4){r) 
with e — [eA + eB)/2 and 5e = (e^ — e_B)/2 where the 
dielectric constant of A (B) compound is denoted as €a 

Macrophase separation is governed by the following set 
of equations for 4){r) and the velocity field v{r). The local 
volume fraction (p obeys 



(4) 



where L is a mobility coefficient. We introduce the po- 
tential U as E{r) = VJ7. The functional derivative of F2 
with respect to U 
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gives us the Maxwell equation 

V{eE) = VpieVpU) = 



(6) 



(the repeated indices iniply summation throughout the 
paper.) The condition 



V X E = 

is automatically satisfied. 

The local velocity field is governed by 



(7) 



p-g^+ p{v-'^)v^-S/p-(l)\7—+'nQ\7\! (8) 

where p is the fluid density, p is determined to satisfy 
the incompressibility condition V • w = 0, and the vis- 
cosity 770 is assumed to be constant and independent 
of </>(r). Equation ([8]) may contain another term like 
UV{SF2/SU). However, this term vanishes because of 
the relation ([5]). 

It is well known that —(t>V{5Fi/5(t)) is related with the 
stress tensor (f) arising from the gradient term in ([2]) 
as 



(9) 



The off-diagonal parts of (r ) are given by 

a^pir) = -X(V„0(r))(V/3(/)(f)) (10) 

In parallel, we may write the electric contribution from 
F2 as 



(11) 



where tT^^(r) is identified with the local Maxwell stress 
defined by 



(12) 



To prove this, let us rewrite the left-hand side of eq. (fTTj) 

as 

f5F2\ , SF2 



Stt 



(V„e)(V/5C/)(V^[/) 



e{V pU){V pU) (13) 



where we have absorbed a term in the form of V X 
into the term Vp in eq. ^ to reach the final expres- 
sion. Then, the relation (|lip follows by calculating the 
right-hand side with the definition of the Maxwell stress 
(eq. m) 



1 



^P<p = T-V^[e(V„C/)(V;3C/)] 



An 
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e{VpU){V^VpU), 



(14) 



where we have used the relation ([6]). 

Macroscopic stress — In principle, the dynamics and 
rheology of the system can be studied by the above set of 
basic equations. However, this usually requires rather in- 
tense numerical computations, and a more coarse grained 
description is called for both to develop the analyti- 
cally tractable theory and to get a deeper physical in- 
sight. In this direction, Doi and Ohta proposed a semi- 
phenomenological rheological constitutive equation for 
the immiscible blends under flow where it was im- 
portant to realize the stress expression in terms of the 
domain configurations, i.e., the stress cr^^ = (o'^^(r)) 
ari sing from the gradient term (eq. (|10p ) can be written 

as [ia [il 



(15) 



Here the bracket indicates the averaging over the spatial 
configurations, F ~ i^/^ is the interfacial tension with ^ 
being the interfacial thickness and the interfacial tensor 
is defined as 



Qaf. 



1 
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dS (rianp) 



(16) 



where J dS () is the surface integral, n is the unit vector 
normal to the interface. In the electro-rheological prob- 
lems, it would be therefore desirable to transform the 
Maxwell stress by treating the F2 term much in the same 
way as the gradient term in the free energy Fi . Below, we 
shall seek such a meaningful expression for the Maxwell 
stress in terms of the external (but not local) electric field 
as well as the domain configurations. 

Equations ([6]) and ^ are solved by a perturbation 
expansion in terms of 5e = (e^ — eB)/2 [l^. The electric 
field is expanded as 



E = + SeE'-^^ + 0(((5e2)) 



(17) 



The external field E'^^ constant in space satisfies the ze- 
roth order solution, i.e., i?*-*^^ = E'^^ . The first order 
solution should satisfy 



v((/)i!;^"=) + eV£;(i) = 

V X = 



(18) 
(19) 



If we put i?*-^^ in terms of the Fourier transformation 
1) 



E^^l = f df Ei^\r)e'''-^ 



5eE^^\ = -—E';''^GZP<Pr. 

P.k g " k 

this satisfies eqs. (fT8|) and ((T9)) . where 



al3 



k - ~gr 



(20) 



(21) 



and its inverse Fourier transformation is given by 

G-^'^f,?) = G(f,P)V"V"^ (22) 
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with -V^G(f, P) = 5{r~r'). Using eq. (HO]), the MaxweU 
stress eq. ([T^ is written up to the order of ((5e)^ as 



Mo 



dr'G'^\r,?)(t>{r)(l){r') 
d?G°'\r,?)4>{r)<t){?) 
dr' I d^'G^^if,?) 



X Gf^^{r,r'')(l){r')(l){r'') 



(23) 



where K, = (,5e)V(47re) and ct^;^' = (e/47r)£;;^^£:|^ is the 
trivial average term, and will be omitted in what follows. 
After averaging over the system volume V , a macro- 
<b) = Iv drM^ is obtained 



scopic stress 
as 



M 
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G^Ji 



^-E^'^Er I GVG^t 



-k 



r {k^Ef + kpEl-){k ■ E'' 
V h 



fc2 



+ 



V 



K, [ kakp{k ■ E<'='f 



k^ 



(24) 



Domain structures — The system governed by eqs. ^ 
and ([S]) generally undergoes macrophase separation in 
which domain coarsening proceeds. There are typically 
two procedures in order to examine a response from 
steady domain structures. One is to study the late stage 
of the phase separation process where droplets make a 
quick response to the electric field compared with the 
coarsening dynamics. The other is to apply steady shear 
flow so that the size of droplets remains finite due to 
their break-up and reconnection. We shall show below 
that in such situations, the Maxwell stress can be evalu- 
ated in terms of the domain structure characteristics. To 
this end, we apply a pre-averaging approximation to the 
double kernel term in eqs. and ([M)) as following: 



f k^kp{k-E^^f , , 



K 



+ {a^ (3) 

^ f {k^Ef + kpE%-){k-E'-) 
2V.h k^ 



where the pre-averaging indicates {kakp/k'^)p,a = 
1/3 6ap and the constant C = 3 is determined by re- 
quiring that the result becomes consistent after taking 



the trace. Substituting eq. f25|) into eqs. (1231) or (|24p . we 
obtain 

^ap Eg 

X / drGf"' (fi ,ri+r) {cj){f, )(/)(fi + r")) 



-t- (a ^ /3) 

= -^^E^J^Er I dfG(f)5f'(f) + (aO/3)(26) 



where we have introduced the pair correlation of - 
V0 as 

g{r) = ((/)(ri)0(fi +r)) 



a3 { ->\ 

9i yr) 



(V?0(fi)vf</>(r2))|,^,=,. 



' and 

(27) 
(28) 



Equation (j26p can be further transformed with the aid 
of the Porod law, which is generally valid in systems 
with domain structures where the interface boundaries 
are rather sharp; R ^ ^ with R being the typical do- 
main scale and the order parameter (t>(f) takes the stable 
value (normalized as 00 = 1) except for the interface re- 
gions [13 • These conditions are satisfied at the late stage 
of phase separation. 

In the length scale £, r ^ R, the structure fac- 
tor V J df g{r)e'^^'^ = {4>^(t>_f:) exhibits the Porod tail 

(0fe0_fc) — Qk~'^, where Q = qaa = J dS/V is the inter- 
face area density. This can be written in the real space as 
51"" (r) ~ 2Q/r where we explicitly write the coefficient 
for the three-dimensional case [l6j . 

Performing the integration in eq. (j26p . we arrive at our 
main result 



(a ^ /3) 

-T{qa5 ssp + qps S5a) 



2q5p 



(29) 



where we cut off the integral at the domain scale, and the 
contribution from the lower bound in the thin interface 
limit ^/i? — 7> is irrelevant. The similarity of cr^jg with 
the usual domain contribution cr^^ (eq. (jlSp) is striking. 
The dimensionless coupling tensor 



ATiK.El^EfR 
f 



(30) 



measures the relative importance of the electrostatic con- 
tribution to the interfacial tension at the scale of the do- 
main size. We see that the Maxwell stress contribution 
leads to an anisotropic renormalization of the interfacial 
tension. 

Comparison with experiment — We now attempt a 
quantitative comparison of our prediction with an exper- 
imental result. Recently, Orihara et. al. have examined 
in detail the behaviors of an immiscible blend subjected 
to a step AC electric field E"^^ = (0, 0, E^^) under steady 
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shear flow Vfiow = (7^,0,0) [2^. In their experiment, 
two polymers with equal viscosity were blended, the mi- 
nority phase of which (with the average volume fraction 
{4'b) = 1/9) forms dispersed droplets in the absence of 
an electric field. Upon the application of the step elec- 
tric field, they observed the droplets elongation along the 
electric field, their coalescence, leading to the network 
structure formation. During this transient process, they 
measured the shear stress axz, and at the same time, 
acquired the three-dimensional images with a confocal 
scanning laser microscope. This allowed them to quanti- 
tatively estimate the interfacial tensor qap and the area 
density Q. In addition, by intermittently turning off and 
on the electric field and the shear flow, they decomposed 
the total shear stress as axz = crY^+a^^+a^^. From these 
measurements, they confirmed the relation = rjoj for 
the viscous stress, a^^ — —Tq^z for the interfacial stress 
(eq. (Uni)) as expected. Furthermore, they found the pro- 
portionality relation between the electric component a^f 
and the interfacial tensor q^z] 



Q 



(31) 



with the slope A^^^p) = 64 (Pa). 

Our formula cq. (|29|) applied to their experimental sit- 
uation reads 



^Qxz'^z 



qxz 



(32) 



where we have used the relation Q ~ {4>b)/R- Substi- 
tuting the experimental parameters {(f) b) = 1/9, E'^^ = 
6/V2 kV/mm (effective value of the AC field), eA = 
2.7eo, CB = 16eo (with eo being the vacuum permit- 
tivity) [l^l, we obtain the theoretical value A'^*''^ = 
4TTKj{E'"^f{(j)B) = 83.7 (Pa) in good agreement with the 
experimental value. This validates the present formula 
eq. (j29p up to a numerical constant, which will enable us 
to investigate rheological response of domains in a sys- 
tematic manner by solving the time-evolution equations 
for the interfacial tensor. 

Summary — It is well known in electrostatics that a 
Maxwell stress is created at the boundary with dielectric 
gap. It would be therefore natural to expect that the 
Maxwell stress in immiscible blends should be correlated 



with the interface configuration. We have demonstrated 
that this is indeed the case, and our formula (eq. (HH)) 
identifies the Maxwell stress contribution as an excess 
interfacial tension. Importantly, this renormalization of 
the interfacial tension is anisotropic, i.e., active only for 
stress components connected with the external electric 
field through the coupling tensor (eq. (PH)) ). 

One may question the accuracy of the pre-averaging 
approximation (eq. ((25|) ). Indeed, the pre-averaging was 
performed in an isotropic state, so we expect that it re- 
tains its physical justification when the magnitude of the 
coupling constant (eq. (PH)) ) is small. Such a situation is 
realized when a strong shear flow is applied in which a 
factor controlling the domain size is the balance between 
surface and viscous stresses 0, i.e., R ~ r/(77o7), where 
7 is the shear rate. Therefore, we find an alternative 
expression for the coupling tensor 



Sa0 



A-KK^El^Ef 



(33) 



This expression implies that the cross coupling between 
the flow and the electric field is highly nontrivial. The 
inverse of its trace ~ "qoj /[4:TTKf{E'^'^)^] is known as 
a Mason number in the electro-rheological literature. 

A remark is in order. The Maxwell stress plays a 
central role in the rheology of magneto-responsive fluids 
where magnetic colloids constitute anisotropic clusters 
such as a chain-like structure under magnetic fleld plj . 
However, to our knowledge, there are no theories to con- 
nect the Maxwell stress with the interfacial tensor as has 
been formulated in the present paper. 

Finally, although semi-quantitative agreement with 
the experiment is encouraging, we have to note that the 
formula (j29p has been derived by assuming that 6e is 
small. With this limitation in mind, we expect that the 
proposed Maxwell stress formula in terms of the interfa- 
cial tensor provides a natural route to construct a coarse 
grained description of the electro-rheology of immiscible 
blends. 
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